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I.  Introduction 

Thermal  fluctuations  and  other  noise  sources  play  an  im¬ 
portant  role  in  the  dynamics  of  conventional  auto-oscillators. 
Thermal  fluctuations  (or  noise)  determine  one  of  the  most 
important,  from  the  practical  point  of  view,  parameters  of  the 
auto-oscillator  -  generation  linewidth  Au>.  The  influence  of 
thermal  noise  is  even  more  important  for  the  nano-scale  spin- 
torque  oscillators.  Due  to  the  extremely  small  sizes  of  STO, 
the  characteristic  operational  energy  of  these  oscillators  can 
be  comparable  with  the  thermal  energy  k-gT  (where  &b  is  the 
Boltzmann  constant  and  T  is  the  thermodynamic  temperature). 
As  a  result,  not  only  the  generation  linewidth,  but  all  the  other 
characteristics  of  spin-torque  oscillators  can  be  dependent  on 
the  temperature.  In  particular,  thermal  fluctuations  lead  to  the 
blurring  of  the  generation  threshold  in  STO  and  to  a  finite 
level  of  power  fluctuations  in  the  above-threshold  operational 
regime,  in  which  the  influence  of  thermal  fluctuations  is, 
usually,  negligible  for  conventional  macro-sized  oscillators. 

In  this  section  we  describe  how  the  influence  of  noise  can 
be  accounted  for  in  the  nonlinear  oscillator  model  (Sec.  II) 
and  consider  several  important  noise-related  issues,  namely, 
the  influence  of  noise  on  the  generation  power  of  an  oscillator 
(Sec.  Ill)  and  the  broadening  of  the  oscillator  linewidth  under 
the  action  of  thermal  fluctuations  (Sec.  IV). 

II.  Stochastic  Nonlinear  Oscillator  Model 

The  stochastic  dynamics  of  an  auto-oscillator  can  be  de¬ 
scribed  in  the  framework  of  the  nonlinear  oscillator  model  (1.5) 
with  external  force  /n(i)  representing  the  action  of  thermal 
fluctuations: 

rj n 

-Jt  +  M\c\2)c  +  r+(|c|2)c  -  T_(|c|2)c  =  /„(f)  .  (1) 

For  any  particular  oscillator  scheme  (1)  can  be  derived  from 
the  stochastic  ’’material”  equations.  For  instance,  to  obtain 
(1)  for  the  Van  der  Pol  oscillator  (see  Sec.  I.IV-A)  one 
has  to  add  to  the  electrical  scheme  shown  in  Fig.  I.2a  an 
additional  stochastic  voltage  source  that  describes  thermal 
Johnson-Nyquist  noise  connected  with  the  resistor  R.  For  a 
spin-torque  oscillator  stochastic  equation  (1)  can  be  derived 
from  the  Landau-Lifshits-Gilbert-Slonczewski  equation  (1.1) 
with  an  additional  random  thermal  magnetic  field. 

For  an  auto-oscillator  of  an  arbitrary  nature,  stochastic  force 
/n(<)  in  (1)  can  be  considered  as  a  phenomenological  term 
describing  action  of  thermal  fluctuations.  In  such  a  case, 
the  statistical  properties  of  fn(t)  should  be  determined  from 
the  condition  of  a  proper  thermodynamical  behavior  of  the 
oscillator  in  the  state  of  thermal  equilibrium  (i.e.,  when  the 
negative  damping  is  vanishing  F_ (|c|2)  =  0).  This  is  achieved 


by  selecting  /n(i)  to  be  a  white  Gaussian  noise  with  the  zero 
mean  and  second-order  correlator  given  by  [1] 

</n(t)/„(t')>  =  0  ,  =  2 DnS(t  -  t')  .  (2) 

Here  Dn  is  the  effective  ’’diffusion  coefficient”  that  character¬ 
izes  the  noise  amplitude.  To  describe  correctly  the  stochastic 
dynamics  of  a  nonlinear  oscillator  with  arbitrary  dependences 
of  the  frequency  u>(p)  and  natural  damping  r+(p)  on  the 
oscillation  power  p,  one  has  to  assume  that  the  diffusion 
coefficient  Dn  also  depends  on  p  as  [1] 

Dn{p)  =  r +(p)v(p)  =  r+(p)^^y  O) 

where  r](p)  =  k^T  /  \uo(jp)  is  the  effective  noise  power  in  the 
nonlinear  regime.  As  it  will  be  shown  below,  this  form  of  the 
diffusion  coefficient  provides  correct  equilibrium  properties 
(see  (7))  for  the  arbitrary  power  dependences  cu(p)  and  r+  (jp) . 

In  (3)  factor  A  is  a  scale  factor  that  relates  the  oscillator 
energy  £{p)  to  the  dimensionless  oscillator  power  p  =  |c|2: 

£{p)  =  A  [  ui{p')dp'  .  (4) 

Jo 

Clearly,  A  depends  on  the  normalization  of  the  oscillator 
amplitude  c.  For  our  choice  of  normalization  of  the  amplitude 
c  of  the  spin-torque  oscillator  (see  (1.15))  A  has  the  form 

A  =  VeffMo/'J 

where  Ves  is  the  effective  volume  of  the  magnetic  material 
of  the  ’’free”  layer  involved  in  the  auto-oscillation,  Mq  is 
the  saturation  magnetization  of  the  ’’free”  layer,  and  7  is  the 
gyromagnetic  ratio. 

For  many  problems,  instead  of  solving  the  nonlinear 
stochastic  Langevin  equation  (1)  for  the  complex  oscillation 
amplitude  c(f),  it  is  convenient  to  consider  the  dynamics  of 
the  probability  density  function  (PDF)  V(t,p,(f>)  that  gives 
the  probability  of  the  auto-oscillator  to  have  the  power  p  and 
phase  (J)  at  the  time  t  (given  some  probability  distribution  at  the 
initial  moment  of  time).  While  the  dynamics  of  the  complex 
amplitude  c(t)  is  described  by  the  stochastic  and  nonlinear 
equation  (1),  the  evolution  of  the  PDF  V  is  governed  by  the 
deterministic  and  linear  Fokker-Planck  equation  that  can  be 
obtained  from  (1)  using  a  standard  formalism  [2] 


dV 

dt 


|[2„(r+-r_)p] 


—  LO 


dV 

d(j> 


d_ 

dp 


d_  (DndV\ 
d(j)  \  2p  d(f> ) 


(5) 


Terms  in  the  left-hand-side  part  of  (5)  originate  from  the 
deterministic  terms  in  (1),  whereas  the  terms  in  the  right-hand- 
side  part  of  (5)  describe  the  stochastic  action  of  the  thermal 
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noise  fn(t).  One  can  see,  that  the  noise-related  terms  in  (5) 
have  the  form  of  a  usual  diffusion  operator  (in  the  power-phase 
coordinates)  with  the  power-dependent  diffusion  coefficient 
Dn.  Thus,  thermal  noise  leads  to  the  diffusion  of  the  oscillator 
in  the  phase  space. 

Equation  (1)  or  (5)  can  be  applied  for  the  analysis  of  a 
number  of  noise-related  problems  of  oscillator  dynamics.  It 
is  clear,  that  one  can  easily  add  to  this  stochastic  scheme 
additional  action  of  deterministic  external  signals  [3]  or  mutual 
coupling  between  oscillators.  It  is,  also,  possible  to  generalize 
(1)  to  include  effects  of  other  sources  of  noise,  e.g.,  in  the 
case  of  spin-torque  oscillators,  flicker  noise  related  to  the 
bias  current  or  mag-noise  related  to  the  fluctuations  in  the 
orientation  of  the  ’’fixed”  layer  magnetization. 

III.  Power  Generated  by  an  Oscillator  in  the 
Presence  of  Noise 

In  this  section  we  consider  the  problem  of  determination  of 
the  power,  generated  by  an  auto-oscillator  in  the  presence  of 
thermal  fluctuations.  In  the  framework  of  the  nonlinear  auto¬ 
oscillator  model  one  can  find  a  complete  analytical  solution 
of  this  problem,  valid  both  below  the  generation  threshold 
(when,  like  in  a  passive  oscillator,  the  oscillations  exist  only 
due  to  the  thermal  noise)  and  in  the  above-threshold  generation 
regime  (when  the  thermal  noise  leads  to  the  fluctuations  of 
the  generated  power).  The  obtained  general  results,  applied 
to  the  case  of  an  STO,  suggest  a  precise  method  for  the 
determination  of  the  threshold  current  It  h  and  quantitatively 
describe  experimental  data  in  the  near-threshold  region,  when 
the  influence  of  the  thermal  fluctuations  can  not  be  ignored. 


A.  Theoretical  Results 


To  find  the  average  power,  generated  by  an  oscillator  in  the 
presence  of  thermal  noise,  it  is  sufficient  to  find  the  stationary 
solution  Vo  of  the  Fokker-Planck  equation  (5).  It  should  be 
noted,  that,  since  (5)  does  not  depend  on  the  oscillator  phase 
<j>  explicitly  (which  is  connected  with  the  phase-invariance  of 
the  Langevin  equation  (1)  or,  in  more  general  terms,  with  an 
autonomous  character  of  the  auto-oscillator’s  dynamics),  the 
stationary  PDF  Vo  is  independent  of  the  phase  <fr.  Thus,  in 
the  stationary  state  any  value  of  the  oscillator  phase  has  equal 
probability,  and  stationary  PDF  Vo  is  a  function  of  only  the 
oscillator  power  p:  Vo  =  Vo  (p)  ■  In  this  case  (5)  becomes  an 
ordinary  second-order  differential  equation 


d 

dp 


2p  (r+  —  r_)  Vo  +  2pDn 


dv 0 

dp 


=  0 


and  a  physically-consistent  solution  of  this  equation  has  the 
form  [1] 


Vo  (p)  =  No  exp 


A 


(6) 

Here  No  is  the  normalization  constant  determined  from  the 
normalization  condition 

noo 


roo 

/  Vo  (p)dp=  1  . 
Jo 


One  can  see  from  (6),  that  at  the  thermal  equilibrium 
(T_(p)  =  0)  stationary  PDF  Voip)  takes  the  form 


Veqip)  =  N0  exp 


A 

hif 


u(p')dp' 


=  No  exp 


£(pY 

kBT_ 


(7) 


which  coincides  with  the  standard  Boltzmann  distribution.  We 
would  like  to  stress,  that  our  approach  gives  a  correct  equi¬ 
librium  distribution  (7)  and,  hence,  satisfies  the  fluctuation- 
dissipation  theorem  [4]  for  an  arbitrary  dependence  of  the 
oscillator’s  frequency  cj(p)  and  damping  T+(p)  on  the  power 
p ,  even  if  these  functions  significantly  change  on  the  thermal 
energy  scale  kBT. 

Equation  (6)  allows  one  to  determine  any  stationary  charac¬ 
teristic  of  a  nonlinear  auto-oscillator  in  the  presence  of  noise. 
For  instance,  the  mean  oscillation  power  p  is  found  to  be 

/>oo 

P=  pVo(p)dp  . 
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The  power  fluctuations  A p2  can  be  calculated  as 


A  p2  =  p2  —  (p)2 


where 

/»oo 

P2  =  /  P2rPo(p)dp  ■ 

Jo 

When  the  temperature  T  is  small,  the  argument  of  the 
exponential  function  in  (6)  is  large  and  (6)  can  be  simplified 
in  two  important  cases.  First,  in  the  below-threshold  regime 
(when  T_(p)  <  T+0)  for  all  p)  the  distribution  Vo(p)  has  a 
maximum  at  p  =  0.  Expanding  the  integrand  in  (6)  in  a  Taylor 
series  near  p  -  0  and  keeping  only  the  first  non-zero  term, 
one  obtains  an  approximate  expression  for  the  PDF  Vo(p)  in 
the  below-threshold  regime: 


V<{p) 


1-C 

r?(0) 


exp 


i-c 

*7(0) 


P 


where  the  supercriticality  is  defined  as  (  =  T_(0)/r+(0), 
which  coincides  with  the  previously  used  notation  (1.26)  for 
the  case  of  a  spin-torque  oscillator.  The  average  power  p  and 
power  fluctuations  A p2  in  this  regime  are  given  approximately 
by  the  expressions: 


P 


Ap2 


*7(0)  =  (kNN\  r+(0) 
1-C  VA^(O);  r+(0)-r_(0) 


(8a) 

(8b) 


The  below-threshold  PDF  V<  (p)  has  approximately  the 
same  Boltzmann-like  form  V<(p)  ~  exp  (—p/p)  as  the  equi¬ 
librium  distribution  Peq(p)  (note,  that  for  small  temperature  T 
one  can  approximate  £{p) /kBT  as  Xco(0)p/kBT  =  p/p{ 0)  in 
(7)),  but  with  the  increased  average  power  level  p.  One  may 
interpret  this  effect  as  an  effective  increase  of  the  temperature 
T  — >  T/(l  —  C)  (see,  e.g.,  [5],  [6]).  We  would  like  to  note, 
however,  that  this  interpretation  is  not  completely  correct,  as 
the  negative  damping  T_(p)  in  this  regime  does  not  only 
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increase  the  power  p,  but,  also,  reduces  the  linewidth  of 
oscillations  (see  Sec.  IV),  which  can  not  be  described  by  any 
changes  of  the  effective  temperature. 

In  the  above-threshold  generation  regime  stationary  PDF 
Vo  ip)  has  a  maximum  at  the  stationary  power  p  =  po  (see 
(1.25)),  for  which  r+(p0)  =  T_(p0),  and  the  integrand  in  (6) 
is  zero.  By  expanding  the  integrand  in  (6)  in  a  Taylor  series 
near  p  =  po,  one  can  derive  an  approximate  expression  for  the 
PDF  Vo  (p)  in  the  above-threshold  regime 


V>{p) 


sff/nAp 


exp 


Here 

A  p2  = 


r+(p0) 


PohiPo) 


(p-Po)2 

2A  p2 

r+(Po)\  kBT  2 


£(Po) 


Po  (9) 


is  the  level  of  power  fluctuations  of  the  oscillator  in  the 
above-threshold  regime  (rp  is  defined  by  (1. 28b)).  The  mean 
oscillation  power  p  in  this  regime  is  approximately  equal  to 
the  zero-temperature  value  po  (1.25). 

One  can  see,  that  the  probability  distribution  function  V>(p) 
in  the  above-threshold  regime  has  a  functional  form  that  is  sub¬ 
stantially  different  from  the  form  of  the  below-threshold  PDF 
V<  (p) .  It  is  also  interesting  to  note,  that  the  power  fluctuations 
in  this  regime  increase  with  temperature  as  A p  ~  y/T,  while 
below  the  generation  threshold  A p  ~  T.  Since  for  typical 
auto-oscillators  Tp  ~  F-(pri),  the  relative  level  of  power 
fluctuations  Ap/po  ~  \A?(po)/po  decreases  with  po  (i.e.,  with 
the  increase  of  the  supercriticality)  and,  for  sufficiently  large 
Po ,  power  fluctuations  become  negligible. 

Neglecting  the  dependence  of  the  oscillator  frequency  u>(p) 
on  the  power  p  (t](p)  =  rj  =  const),  which  has  only  a 
minor  influence  on  the  form  of  the  PDF,  and  using  (1.20), 
one  can  derive  from  (6)  an  explicit  general  expression  for 
the  stationary  probability  distribution  function  of  a  spin-torque 
oscillator  valid  in  the  whole  range  of  bias  currents: 


Vo{p)  = 


exp  [-(1  -  C  +  Cp)2/2(r]\ 
nr]  1  +  erf  [(C  -  l)/VKv\ 


for  Q  =  0  and 


Vo(p)  = 


Qexp  [-«  +  Q)(l  +  Qp)/Q2v] 
(l  +  Qp)PE0(((  +  Q)/Qiv) 


(10a) 


(10b) 


for  Q^O.  Here  /3  =  —(1  +  Q)C,/Q2r],  erf(a;)  = 
{2/yfn)j^e~t  dt  is  the  error  function,  and  En(x)  = 


ft 


t/tndt  is  the  exponential  integral  function.  Fig.  1 


shows  the  profile  of  the  probability  distribution  function  (10) 
in  the  below-threshold  (£)  =  0),  threshold  ((  =  1),  and  above¬ 
threshold  (C  =  2)  regimes. 

Using  the  analytical  expression  (10)  for  the  probability 
distribution  function,  it  is  possible  to  obtain  an  analytical 
expression  for  the  mean  oscillation  power  p  of  a  spin-torque 
oscillator: 


_=  f2p  exp  [— (C  -  1)2/2C??]  C  -  1 

P  V  TfC  1  +  erf  [(C  —  1)/a/2C^7]  C 


(11a) 


Fig.  1.  Stationary  probability  distribution  function  Vo  (p)  (10)  for  several 
values  of  the  supercriticality  parameter  £  (from  [1]).  Nonlinear  damping 
coefficient  Q  =  0.3,  noise  level  77  =  0.05.  Data  for  £  =  0  (thermal 
equilibrium)  have  been  divided  by  2. 


Fig.  2.  Dependence  of  the  mean  oscillator  power  p  (1 1)  on  the  supercriticality 
f  for  several  noise  levels  rj  (from  [1]).  Nonlinear  damping  coefficient  Q  = 
0.3.  The  lower  curve  for  77  =  0  coincides  with  the  deterministic  value  (1.25) 
of  the  oscillator  power  po . 


for  Q  =  0  and 


Qi  1  L  exp(— (C  +  Q)!Q2rt)' 

Q  +  C  L  E0(((  +  Q)/Q2r1)  _ 


c-i 

C  +  Q 


(lib) 


for  Q  ^  0.  Fig.  2  shows  the  dependence  (11)  of  the  mean 
power  p  on  the  supercriticality  £  for  several  values  of  the 
noise  power  rj. 


B.  Comparison  With  Experiments 

To  illustrate  the  application  of  the  above  described  theory  to 
a  real  experimental  situation  we  compare  below  the  theoretical 
expressions  for  the  auto-oscillator  power  in  the  presence  of 
noise  with  the  recent  measurements  of  the  generated  power  in 
a  GMR  spin-valve  auto-oscillator  based  on  a  metallic  nano¬ 
pillar  [7].  The  experimental  data  from  [7]  are  shown  by  black 
dots  in  the  main  frame  of  Fig.  3.  To  analyze  these  experimental 
data  it  is  important,  first  of  all,  to  determine  the  threshold  of 
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Fig.  3.  Main  panel:  Dependence  of  the  mean  power  p  on  the  bias  current 
I.  Dots  -  experiment  [7]  for  T  =  225  K,  solid  line  -  theoretical  dependence 
(11)  for  Jtn  =  4.9  mA,  p  =  4.2- 10~4,  Q  =  0.3  (from  [1]).  Inset  shows  the 
same  data  for  inverse  power  1  /p  in  near-threshold  range  of  currents.  Dashed 
line  corresponds  to  the  approximate  expression  (12)  valid  for  small  currents. 
Intersection  of  this  line  with  x  axis  gives  the  value  of  the  threshold  current 
Ah. 

microwave  generation,  i.e.  the  magnitude  of  the  bias  current  at 
which  the  spin-torque  oscillator,  instead  of  passively  filtering 
its  eigen-frequency  from  the  thermal  noise,  starts  to  actively 
self-generate.  We  believe  that  this  can  be  done  by  analyzing 
the  experimental  data  for  very  small  values  of  the  bias  current. 
In  the  inset  of  Fig.  3  we  show  by  dots  the  experimental 
dependence  of  the  inverse  mean  power  1  /p  on  the  bias  current 
I.  These  dots  can  be  fitted  to  the  theoretical  expression  (8a) 
which  gives  a  linear  dependence  of  the  inverse  power  on  the 
bias  current  for  small  values  of  current 

=  oc  (1th  -  1)  ■  (12) 

V 

This  linear  dependence,  taking  place  in  the  below-threshold 
regime  I  <  It h  and  shown  by  a  dashed  line  in  the  inset  of 
Fig.  3,  crosses  the  horizontal  axis  (axis  of  the  bias  current)  at 
the  point  I  =  lt h  and,  therefore,  allows  one  to  determine  the 
lth  with  a  reasonably  high  precision. 

Thus,  the  precise  measurement  of  the  oscillator  output 
power  as  a  function  of  the  bias  current  for  low  values  of  current 
in  combination  with  the  theoretical  formula  (8a)  provides  a 
simple  method  for  the  precise  determination  of  the  threshold 
current  lth  of  microwave  generation  in  a  spin-torque  oscillator 
in  situations  when  the  influence  of  thermal  fluctuations  is 
strong,  and  the  determination  of  lth  by  other  means  is  difficult. 
In  particular,  from  the  inset  of  Fig.  3  one  can  immediately 
determine  the  value  of  the  threshold  current  of  f\t  =  4.9  mA, 
which  is  substantially  larger  than  the  current  1*  «  4.2  mA,  at 
which  the  thermally-induced  oscillations  become  observable 
in  the  experiment  [7]. 

Using  thus  determined  value  of  the  threshold  current  lt h  = 
4.9  mA  to  calculate  the  supercriticality  parameter  £  =  I /lth 
in  the  theoretical  expression  (11),  it  is  possible  to  describe 
theoretically  the  power  p,  generated  by  an  STO,  in  a  wide 
range  of  bias  currents  (see  main  panel  in  Fig.  3).  We  would 
like  to  stress,  that  the  correct  determination  of  the  threshold 


Fig.  4.  Main  panel:  Dependence  of  the  inverse  power  1/p  on  the  bias 
current  I.  Dots  -  experiment  [8]  for  a  magnetic  tunnel  junction  (the  direction 
of  the  bias  current  was  changed  to  coincide  with  our  sign  agreement),  solid 
lines  -  theoretical  dependences  (12)  demonstrating  different  values  of  the 
threshold  current,  determined  from  the  data  for  positive  (/th,  [)0s  1  and  negative 
(7th,  neg)  current  directions.  Inset  shows  the  same  set  of  data  for  the  direct 
power  p. 

current  /th  (and,  therefore,  the  parameter  Q  is  essential  for  the 
comparison  of  (11)  with  experiment,  since  the  dependence  of 
the  oscillator  power  p  on  the  other  parameters  of  the  model 
(like  noise  power  rj  and  nonlinear  damping  parameter  Q)  in 
the  above-threshold  region  (£  >  1)  is  relatively  weak. 

Fig.  4  demonstrates  another  example  of  the  threshold  current 
determination  using  the  linear  dependence  (12)  in  the  below- 
threshold  region.  In  this  case  the  above  described  method  of 
the  threshold  current  determination  is  applied  to  a  spin-torque 
device  based  on  a  magnetic  tunnel  junction  (MTJ).  Dots  in 
the  main  panel  of  Fig.  4  show  the  experimental  data  for  the 
inverse  power  1  /p  measured  in  the  below-threshold  regime  in 
an  MTJ  nano-pillar  (see  [8]),  in  which  the  spacer  between  the 
’’free”  and  ’’fixed”  magnetic  layers  (see  Fig.  I.lb)  was  made 
of  a  non-magnetic  insulator  AI2O3.  Due  to  a  much  larger 
resistance  of  the  MTJ  nano-pillar  (compared  to  fully-metallic 
GMR  spin  valve  nano-pillars,  where  the  spacer  is  made  of  a 
non-magnetic  metal),  the  authors  of  [8]  were  able  to  measure 
the  power  of  the  thermally-induced  oscillations  of  the  device 
even  for  the  negative  values  of  the  bias  currents  I  <  0,  when 
the  spin-transfer  torque  creates  effective  positive  damping  and, 
according  to  (12),  reduces  the  measured  oscillation  power  p. 

One  can  see  from  Fig.  4  that  the  experimental  dependences 
of  the  inverse  output  power  on  the  bias  current  for  positive 
I  >  0  and  negative  I  <  0  current  directions  are  both  linear 
(12),  but  have  substantially  different  slopes.  Respectively, 
values  of  the  threshold  current  Jth,  estimated  using  (12) 
from  the  experimental  data  for  positive  (/th,  pos  =  9  mA) 
and  negative  (lth,  neg  =  2  mA)  currents,  are  significantly 
different.  Note,  that  both  /th,  pos  and  /th,neg  correspond  to 
the  same  magnetization  configuration,  but  were  calculated  for 
different  current  polarities.  Since  the  threshold  current  (1.23)  is 
inversely  proportional  to  the  spin-polarization  efficiency  e,  this 
difference  in  Jth  suggests  that  the  spin-polarization  efficiency 
£  for  MTJ  depends  not  only  on  the  mutual  orientation  of  the 
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magnetizations  of  the  ’’free”  and  "fixed”  layers,  but,  also,  on 
the  direction  of  the  bias  current.  In  the  case  shown  in  Fig.  4  the 
ratio  of  spin-polarization  efficiencies  for  positive  and  negative 
currents  is  ep0s/£neg  =  2/9  «  0.22. 

Thus,  the  measurements  of  the  current  dependence  of  the 
noise-induced  output  power  in  MTJ-based  STO  can  provide 
direct  information  about  the  spin-polarization  efficiency,  and 
can  be  used  for  verification  of  theories  of  spin-dependent 
tunneling  in  magnetic  nano-structures. 

We  would  like  to  mention,  that  the  threshold  current, 
estimated  in  [8]  using  the  measurements  of  the  linewidth 
and  the  data  for  the  direct  (rather  than  inverse )  power  p,  is 
Ith  =  20  ±  5  mA,  which  is  far  above  our  estimates.  We 
attribute  this  difference  to  the  fact  that  other  sources  of  noise 
can  give  a  substantial  contribution  to  the  linewidth  of  a  highly- 
resistive  MTJ-based  device,  and,  therefore,  the  estimation  of 
the  threshold  current  based  on  the  linewidth  data  could  be 
unreliable.  On  the  other  hand,  the  dependence  of  the  direct 
power  p  on  the  bias  current  I  (see  inset  in  Fig.  4)  is  nonlinear, 
and  it  is  much  more  difficult  to  make  an  estimate  of  the 
threshold  current  based  on  such  a  nonlinear  dependence.  Thus, 
the  difference  in  the  slopes  for  currents  of  different  polarity, 
that  is  clearly  seen  in  Fig.  4,  which  demonstrates  the  current 
dependence  of  the  inverse  power,  is  hardly  noticeable  in  the 
inset  of  Fig.  4,  where  the  current  dependence  of  the  direct 
power  is  presented. 

IV.  Generation  Linewidth  of  an  Auto-Oscillator 

Thermal  fluctuations  do  not  only  blur  the  stationary  dis¬ 
tribution  Vo(p)  of  the  oscillation  power,  but,  also,  alter  the 
dynamical  behavior  of  an  auto-oscillator.  In  the  presence  of 
noise  the  auto-oscillation  amplitude  c(t)  is  not  anymore  a 
harmonic  function  of  time  c(t)  ~  corresponding  to 

a  simple  monochromatic  oscillation,  but  is,  instead,  a  rather 
complicated  stochastic  function,  and  its  spectrum 


has,  strictly  speaking,  non-zero  spectral  components,  corre¬ 
sponding  to  every  frequency  Cl. 

The  power  spectrum  of  the  auto-oscillation  S{  12),  that  can 
be  defined  as 

(cqCq/)  =  27 tS(Q  —  Cl')S(Cl) 

is  a  convenient  quantity  that  gives  a  detailed  description  of 
the  auto-oscillation  spectrum  in  the  presence  of  noise.  Here  the 
angular  brackets  (...)  stay  for  the  averaging  over  the  statistics 
of  thermal  fluctuations.  For  a  coherent  monochromatic  oscil¬ 
lation  with  the  frequency  oj„  the  power  spectrum  in  singular 
iS(fi)  ~  6(Cl  —  u)g).  In  the  presence  of  noise  S( Cl)  has  a  finite 
width  (called  sometimes  generation  linewidth)  A u>,  which 
provides  a  simple  way  for  a  single-parameter  characterization 
of  a  noisy  auto-oscillation  spectrum.  The  generation  linewidth 
A u>  is  one  of  the  most  important  parameters  of  an  auto¬ 
oscillator,  especially  with  regards  to  its  practical  applications. 

Alternatively,  the  coherence  of  oscillations  can  be  charac¬ 
terized  by  an  autocorrelation  function 

/C(t)  =  { c(t  +  r)c*(t ))  . 


The  descriptions  of  noisy  oscillations  by  a  power  spectrum 
S(Cl)  or  an  autocorrelation  function  /C(r)  are  fully  equiva¬ 
lent,  since  these  two  functions  are  connected  by  the  Fourier 
transform 

S(Cl)  =  J  K,{r)einTdT  . 


As  it  will  be  shown  below,  in  many  cases  the  power  spec¬ 
trum  of  an  auto-oscillation  can  be  approximately  represented 
in  a  simple  Lorentzian  form 


s{ci) 


2A  up 

(0-w)2  +  Aw2 


(13a) 


where  p  is  the  average  oscillation  power,  and  to  is  the  mean 
oscillation  frequency.  Note,  that  in  (13a)  and  in  the  following 
text  Au>  is  defined  as  /;«//-! inewidth  at  a  half-power  level.  The 
autocorrelation  function,  corresponding  to  (13a),  has  a  simple 
exponential  form 


K,{t)  =  p exg>{—itjjT  —  Alo\t\)  .  (13b) 


It  is  important  to  stress,  that  in  many  cases  the  experimental 
data  for  the  generation  linewidth  are  fitted  to  the  simple 
Lorentzian  function  (13),  even  if  the  actual  lineshape  S(Cl) 
(or  autocorrelation  function  /C(r))  has  more  complicated  form. 
Thus,  usually,  the  generation  linewidth  Ac o  is  understood  in  a 
sense  of  an  experimentally  determined  quantity  obtained  as  a 
best  fitting  parameter  in  the  Lorentzian  expression  (13). 

In  the  next  section  we  derive  expressions  for  generation 
linewidth  A u>  of  an  auto-oscillator  is  several  important  partic¬ 
ular  cases,  using  the  stochastic  Langevin  equation  (1). 


A.  Below-Threshold  Regime 

In  the  below-threshold  regime  (T_(p)  <  T+(p)),  when 
the  characteristic  oscillation  power  p  ~  77  is  small,  one  can 
neglect  all  the  nonlinearities  in  (1).  In  such  a  case  equation 
(1)  becomes  an  equation  describing  a  usual  linear  oscillator  in 
the  presence  of  noise: 

^  +  iu(o)c+  [r+(o)  -r_(o)]c  =  /n(f) . 

One  can  easily  find  a  solution  of  this  linear  equation  in  the 
form 

t 

c(f)  =  J  /„(/')  exp{— [iw(0)  +  r+(0)  -  r_(0)](t  -  t')}dt' 

—  OO 

Using  the  statistical  properties  of  the  thermal  fluctuations  (2), 
it  is  easy  to  show  that  the  autocorrelation  function  /C(r)  in 
this  linear  case  has  a  simple  exponential  form  (13b)  with  the 
average  power  p  given  by  (8a),  mean  frequency  defined  by 
u>  =  w(0),  and  the  full  linewidth  given  by 

2A w  =  2r+(0)  -  2r_(0)  =  2r+(0)(l  -  C)  •  (14) 

It  is  clear,  that  the  power  spectrum  S(CT)  of  the  oscillation  in 
this  linear  case  has  exactly  the  Lorentzian  shape  (13a). 

For  a  spin-torque  oscillator  the  expressions  for  the  positive 
and  negative  damping  are  given  by  (1.20),  and  the  linewidth 
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2A u)  of  such  an  oscillator  decreases  linearly  with  the  bias 
current  I  as 

2Alo  =  2rG  -  2cr I  =  2rG  ^1  -  .  (15) 

The  measurements  of  the  current  dependence  of  the  linewidth 
Aw  in  the  below-threshold  regime  provide  an  alternative 
method  of  the  experimental  determination  of  the  threshold 
current  /ti,  for  current-induced  excitations.  As  it  was  discussed 
above,  this  method  may  not  be  very  accurate  for  highly- 
resistive  STO  based  on  magnetic  tunnel  junctions  (MTJ). 

Recently,  the  linear  dependence  (15)  of  the  linewidth  Aw 
on  the  bias  current  I  has  been  observed  both  in  fully  metallic 
GMR  spin  valves  [9]  and  in  MTJ  based  on  the  TMR  effect 
[8], 

B.  Above-Threshold  Regime:  ” Linear ”  Oscillator 

Far  above  the  generation  threshold  the  power  fluctuations 
in  an  auto-oscillator  are  negligible  compared  to  the  mean 
generated  power.  In  this  regime  the  complex  amplitude  c(t) 
becomes  a  pure  phase-modulated  process  c(f)  «  y/poe1^^. 
Consequently,  the  generation  linewidth  Aw  and  the  overall 
shape  of  the  auto-oscillation  power  spectrum  S(Cl)  is  de¬ 
termined  only  by  the  phase  fluctuations.  This  property  is 
general  for  auto-oscillators  of  any  nature,  and  is  connected 
with  different  types  of  stability  of  the  oscillator  dynamics  with 
respect  to  the  power  and  phase  fluctuations  (see  Sec.  I.V-A2). 

Thus,  to  find  a  generation  linewidth  of  an  auto-oscillator 
in  the  above-threshold  generation  regime,  it  is  sufficient  to 
determine  the  statistical  properties  of  the  auto-oscillator  phase 
(pit).  Here  we  will  consider  this  problem  for  a  ’’linear”  oscil¬ 
lator,  i.e.  for  an  auto-oscillator  with  frequency  independent 
of  the  oscillation  power  ui(jp)  =  cug  =  ojq  =  const.  As 
it  was  noted  before,  the  majority  of  conventional  oscillators 
belong  to  this  class,  whereas  spin-torque  oscillators,  in  general, 
are  characterized  by  a  strong  dependence  of  the  oscillation 
frequency  oj(p)  on  the  power  p  (this  case  will  be  considered 
in  the  following  section). 

Using  (1),  one  can  derive  a  stochastic  equation  for  the  phase 
<p  =  arg(c)  of  an  auto-oscillator 

^+Ut=-LlmlUt)]  (16) 

where  fn(t)  =  /n(t)e-l^W.  Note,  that  the  statistical  properties 
of  the  function  fn(t)  are  the  same  as  for  the  function  fn(t),  i.e. 
fn(t)  is  a  white  Gaussian  stochastic  process  with  the  second- 
order  correlator  given  by  (2). 

Far  above  the  generation  threshold  one  can  substitute  the 
stationary  power  po  for  p  in  the  right-hand-side  part  of  (16), 
thus  reducing  this  equation  to  a  closed-form  equation  for  the 
phase  0,  which  has  the  following  solution  for  the  phase  as  a 
function  of  time: 

(p{t)  =  0(0)  -  u)gt  +  — L  [  Im [fn(t')]dt'  , 

\/P o  Jo 

where  it  was  assumed  that  the  phase  at  the  initial  moment 
of  time  t  =  0  has  fixed  value  0(0).  This  solution  describes 
the  ’’Brownian  motion”  of  the  phase  0(f)  under  the  action  of 


thermal  fluctuations.  It  is  clear,  that  0(f)  is  a  Gaussian  process 
with  the  mean  value 


(0(f))  =  0(0)  -Wgf  (17a) 

and  the  variance 

A02(f)  =  (02(f))  -  [(0(f))]2  =  r+(P°MPo)  |f |  •  (17b) 

Po 

The  phase  variance  A 02(f)  increases  linearly  with  the  time 
interval  |f|,  exactly  as  the  square  of  a  particle  displacement 
for  a  usual  Brownian  motion. 

The  autocorrelation  function  /C(f),  corresponding  to  the 
Gaussian  random  phase  0(f),  can  be  written  as 


£{t)  =  Po  (exp{i[0(f)  -  0(0)]}) 

=  exp[—  A0(f)2/2]  . 


For  a  ’’linear”  auto-oscillator  (A02(f)  ~  |f|,  see  (17b))  the 
auto-correlation  fuction  /C(f)  has  the  Lorentzian  form  (13b) 
with  the  average  frequency  u  =  and  full  linewidth 


2Awo  — 


A02(f) 

1*1 


r+(po) 


vjpo) 

Po 


The  ratio  of  the  noise  power  rj(po)  to  the  stationary  auto¬ 
oscillation  power  po  in  the  above  derived  linewidth  expres¬ 
sion  can  be  rewritten  as  a  ratio  of  corresponding  ener¬ 
gies  kpJT / £(po),  yielding  a  physically  transparent  expression 
for  the  Lorentzian  generation  linewidth  of  a  ’’linear”  auto¬ 
oscillator: 


2Acuq  — 


r+(po) 


kBT 

£{po) 


(18) 


To  the  best  of  our  knowledge,  the  generation  linewidth  in  a 
general  form  (18)  was  first  written  in  [10]. 

The  damping  rate  r+(po)  in  (18)  determines  the  over¬ 
all  scale  of  the  possible  linewidth  variations,  and  the  ratio 
kBT  /  £{po)  describes  the  linewidth  reduction  due  to  a  smaller 
influence  of  noise  (having  the  constant  energy  kBT )  for  above¬ 
threshold  oscillations,  the  energy  £(po)  of  which  increases 
with  the  oscillation  power  po .  Note,  also,  that  if  we  formally 
use  (18)  for  the  evaluation  of  the  linewidth  in  the  state  of  ther¬ 
modynamic  equilibrium  (£(po)  — >  kBT,  F+(po)  — ►  r+(0)), 
we  obtain  the  linewidth  value  that  is  2  times  smaller  than  the 
value  given  by  (14)  in  the  same  limit  (r_(0)  — >  0).  This 
difference  is  explained  by  the  fact  that  only  phase  fluctuations 
lead  to  the  linewidth  broadening  of  an  auto-oscillator  in  the 
above-threshold  generation  regime,  whereas,  for  a  passive 
oscillator,  both  phase  and  amplitude  fluctuations  play  the  same 
role. 

The  expressions  for  the  generation  linewidth  of  many  dif¬ 
ferent  types  of  conventional  (’’linear”)  auto-oscillators  have 
been  derived  previously  (see,  e.g.,  a  classical  book  [11]).  It  is 
important  to  note,  that  all  these  expressions  can  be  rewritten 
in  the  general  form  (18).  In  particular,  the  full  generation 
linewidth  of  an  electrical  auto-oscillator  (similar  to  the  one 
considered  in  Sec.  I.IV-A)  is  given  by  the  classical  formula 
[11]: 


2Awo  — 


kBTR  Wg 

r/2 

uo 
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where  R  is  the  resistance  of  the  circuit,  u>o  =  1/y/LC  is  the 
resonance  frequency  of  the  electrical  circuit  having  inductance 
L  and  capacitance  C,  and  L'(l  is  the  voltage  amplitude  of 
self-sustained  electrical  oscillations  on  the  capacitor  C.  This 
expression  can  be  exactly  rewritten  in  the  form  (18),  taking 
into  account  that,  for  an  electrical  oscillator,  the  damping  rate 
is  r+  =  R/(2L ),  and  the  oscillation  energy  is  £  =  CUq/2. 

The  general  expression  (18)  for  the  oscillator  linewidth, 
sometimes,  can  be  used  even  for  the  description  of  quantum 
systems.  For  example,  the  quantum  limit  for  the  linewidth 
of  a  single-mode  laser  (Schawlow-Townes  limit)  leads  to  the 
following  expression  [12],  [13] 


2Aw0  — 


huj0  rg 

2  Pout 


Here  ojq  is  the  generation  frequency,  r0  is  the  loss  rate,  which 
is  assumed  to  be  caused  only  by  the  energy  radiation  out  of 
the  laser  cavity  (losses  in  the  active  medium  were  neglected), 
and  Pout  is  the  output  power  of  the  laser.  This  expression, 
also,  can  be  rewritten  in  the  form  (18),  if  one  uses  the  energy 
of  quantum  fluctuations  Tiloq/2  instead  of  the  thermal  noise 
energy  k\{l\  and  takes  into  account,  that  the  output  power 
Pout  of  the  laser  is  connected  with  the  oscillation  energy  £ 
by  the  expression  Pout  =  ToP 

An  attempt  to  calculate  the  generation  linewidth  of  a  spin- 
torque  oscillator  has  been  undertaken  in  [14],  and  resulted 
in  an  expression  (see  equation  (28)  in  [14]),  that,  also,  can 
be  rewritten  in  the  form  (18).  However,  the  dependence  of 
the  generated  frequency  of  the  generation  power  was  ignored 
in  [14],  and,  therefore,  the  linewidth  expression  obtained  in 
[14]  is  correct  only  qualitatively.  Thus,  it  is  necessary  to 
consider  the  case  of  a  substantially  "nonlinear”  auto-oscillator 
(i.e.  an  oscillator  with  a  large  nonlinear  frequency  shift)  to 
quantitatively  describe  the  generation  linewidth  of  an  STO. 


C.  Above-Threshold  Regime:  ’’Nonlinear”  Oscillator 

The  expression  for  the  auto-oscillator  generation  linewidth 
(18),  which  was  verified  for  many  physical  realizations  of 
auto-oscillators,  underestimates  the  generation  linewidth  of  an 
STO  by  one-two  orders  of  magnitude.  The  reason  for  such  a 
large  difference  is  that  in  the  derivation  of  (18)  the  oscillation 
frequency  was  assumed  to  be  independent  of  the  oscillation 
power  uj(p)  =  const.  For  spin-torque  oscillators,  as  it  was 
already  emphasized,  the  influence  of  the  nonlinear  frequency 
shift  is  very  important.  To  find  the  correct  expression  for  the 
generation  linewidth  of  an  STO,  as  well  as  for  any  other 
"nonlinear”  oscillator,  one  has  to  explicitly  take  into  account 
the  power  dependence  of  the  frequency  u>  (p)  in  the  theory. 

In  the  above-threshold  generation  regime  the  power  fluc¬ 
tuations  8p(t)  =  p(t)  —  po  are  small  even  for  a  ’’nonlinear” 
oscillator  6p(t)  <C  po-  However,  to  obtain  a  correct  result  for 
the  generation  linewidth  one  should  not  ignore  them  altogether, 
but,  instead,  should  use  this  smallness  to  linearize  (1)  near  the 
stationary  generation  power  p  =  po  to  derive  an  approximate 
system  of  equations  for  the  power  fluctuations  6p(t)  and  the 


phase  </(£)  in  a  ’’nonlinear”  auto-oscillator: 

C^+2Tv5p  =  2v/p0Re[/„(f)]  (19a) 

-J7  +U g  =  -7=Im[/n(f)]  -  NSp  (19b) 

at  v Po 

where,  as  in  the  previous  subsection,  /„(£)  =  fn(t)e~ is  a 
stochastic  process  with  the  same  properties  as  /n(f).  The  left- 
hand  side  of  (19a)  is  the  same  as  (1.27)  for  power  fluctuations 
Sp  of  free  auto-oscillations. 

One  can  see  from  (19b),  that  the  frequency  nonlinearity  N 
creates  an  addition  noise  term  — NSp  in  the  the  phase  equation 
(19b).  The  physical  mechanism  behind  this  additional  noise  is 
clear:  for  a  nonlinear  auto-oscillator  the  power  fluctuations 
Sp(t)  lead  to  the  frequency  modulation  uj(p(t))  =  u>(po  + 
5p(t))  ~  Lug  +  NSp.  In  some  sense  this  additional  noise  term 
can  be  considered  as  an  inhomogeneous  broadening  of  the 
oscillator  linewidth:  oscillators  with  different  (due  to  thermal 
fluctuations)  powers  p  have  different  oscillation  frequencies 
uj(p). 

Since  the  stochastic  system  (19)  is  a  linear  system  of  equa¬ 
tions,  and  the  noise  /„(£)  is  a  Gaussian  process,  both  Sp(t)  and 
</>(£)  are  the  Gaussian  processes  also.  One  can  easily  obtain 
a  complete  set  of  statistical  characteristics  of  these  processes. 
Namely,  the  mean  value  of  the  power  fluctuations  following 
from  (19)  is  zero  (Sp(t))  =  0,  whereas  its  correlation  function 
has  the  form 

(■ Sp(t)5p(t '))  =  r+r(Po)  ?KPo)poe~2rpl*~*'1  ■  (20) 

1  p 

The  mean  value  of  the  phase  4>(t)  (for  a  fixed  value  <j>(0)  of 
the  phase  at  the  initial  moment  of  time  t  =  0)  is  given  by 


(</(£))  =  4>(0)  —  Logt  (21a) 

and  the  variance  of  the  phase  fluctuations  is  expressed  as 

m  ,  o  1  —  e_2rpl4l  ’ 


A02(i)  =  2A  w0 


(1 


2rr 


(21b) 


where  v  is  the  dimensionless  nonlinear  frequency  shift  coef¬ 
ficient  (1.34). 

In  contrast  to  the  case  of  a  ’’linear”  oscillator  (17b),  the 
variance  A<f>2(t)  (21b)  has  a  nonlinear  dependence  on  the  time 
interval  |t|.  This  is  caused  by  a  finite  correlation  time  (see 
(20))  of  the  additional  ’’nonlinear”  noise  term  — NSp(t )  in 
(19b).  Only  for  the  time  intervals  \t\,  that  are  much  larger  than 
the  correlation  time  1  /Tp  of  the  power  fluctuations,  the  phase 
variance  A (f>2(t)  starts  to  grow  linearly  with  |£|.  Consequently, 
the  spectrum  of  a  ’’nonlinear”  auto-oscillator  is,  in  general, 
non-Lorentzian.  There  are,  however,  two  limiting  cases  when 
a  rather  complicated  expression  (21b)  for  the  phase  variance 
of  a  nonlinear  auto-oscillator  can  be  substantially  simplified. 

If  the  temperature  T  and,  consequently,  the  generation 
linewidth  A u>  is  sufficiently  small  Ac o  <C  I’p,  one  can  neglect 
the  exponential  factor  in  the  second  term  of  the  right-hand- 
side  part  of  (21b)  at  the  characteristic  decoherence  time  scale 
of  |f|  ~  1/At j.  Then,  the  following  expression  for  the  phase 
variance  is  obtained  in  the  low-temperature  limit: 

A</>2(f)  ~  2Acu0(l  +  v2)\t\  —  . 
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In  this  limit  the  autocorrelation  function  JC(t)  becomes  expo¬ 
nential  (13b),  and  the  power  spectrum  S(Q)  has  a  Lorentzian 
shape  (13a)  with  the  full  linewidth  of  generation  given  by  the 
expression 


2 Acu  —  2Acuq(1  T  zz“) 


(i  +  v2)r+{Po) 


kBT 

£{po) 


(22) 


One  can  see,  that  the  nonlinearity  of  the  auto-oscillator  fre¬ 
quency  uj(p')  leads  to  the  increase  of  the  generation  linewidth 
by  (1  +i/2)  times.  The  expression  for  the  generation  linewidth 
of  a  ’’nonlinear”  auto-oscillator  (22)  was  derived  on  the  basis 
of  a  nonlinear  oscillator  model  and  published  for  the  first  time 
in  [10], 

One  can  show,  that  the  condition  Ac o  Tp  of  applicability 

of  (22)  can  be  rewritten  as  a  condition  for  the  temperature 


kBT  <C 


£Qo) 
1  +  v2 


(23a) 


or,  using  (9),  as  a  condition  for  the  relative  power  of  thermal 
fluctuations 


A  p 

—  < 
Po 


Vl 


(23b) 


Thus,  (22)  represents  a  low-temperature  limit  of  the  gen¬ 
eration  linewidth  of  a  ’’nonlinear”  auto-oscillator.  The  esti¬ 
mations  for  the  typical  parameters  of  a  spin-torque  oscillator 
(permalloy  circular  nano-pillar  of  the  radius  Rc  =  50  nm 
and  thickness  of  the  ’’free”  layer  L  =  5  nm,  generation 
frequency  cug/27r  ~  30  GHz)  show  that  the  expression  (22)  is 
quantitatively  correct  for  the  temperatures  T  <  10  —  100  K, 
depending  on  the  supercriticality  parameter  £. 

Another  limiting  case,  which  allows  an  easy  analysis  of 
(21b),  is  the  opposite  limiting  case  of  relatively  long  (com¬ 
pared  to  the  inverse  linewidth  1/At u)  correlation  times  of 
the  power  fluctuations  1/Tp,  or  the  case  of  relatively  large 
generation  linewidths  Acu>Tp.  In  this  case  the  exponential 
function  in  (21b)  can  be  developed  in  a  Taylor  series  to  give: 


A cj)2{t)  k,  2Aw0(|f|  +  v2Vpt2)  . 


For  sufficiently  large  frequency  nonlinearities  v  1  one  can 
drop  the  first  (linear  in  |f|)  term  in  the  brackets,  and  retain 
only  the  second  term  that  is  one  quadratic  in  time.  Then,  the 
phase  variance  A 4>2(t)  ~  t 2,  and  the  power  spectrum  takes 
the  form 


S(VL)  ~  exp 


(cu(po)  -  O)2' 

2  Aw*2 


with  the  characteristic  linewidth 


In  this  limit  the  auto-oscillator  linewidth  is  large  Aw* 
rp  and  proportional  to  the  normalized  nonlinear  frequency 
shift  coefficient  v  in  the  first  power.  Also,  in  contrast  with  the 
low-temperature  linewidth  expression  (22),  giving  the  linear 
dependence  on  temperature,  the  expression  (24)  for  the  auto¬ 
oscillator  linewidth  is  proportional  to  \/T.  It  is  interesting  to 
note,  that  formula  for  the  auto-oscillator  linewidth,  which  is 
qualitatively  similar  to  (24),  was  for  the  first  time  proposed 


in  [15]  for  the  case  of  a  spin-torque  auto-oscillator  (see  (2)  in 
[15]).  In  the  same  paper  the  \jT  temperature  dependence  of 
an  STO  linewidth  was  obtained  as  a  result  of  direct  numerical 
simulations. 

The  expression  (24)  for  the  generation  linewidth  is  valid  in 
the  temperature  range 

(r^))^«^«(r^))£(»»>  (25‘> 

or,  in  terms  of  power  fluctuations, 

1  A» 

i-r  <  —  <  1  •  (25b) 

M  Po 

Thus,  (24)  is  the  asymptotic  expression  of  the  generation 
linewidth  of  a  strongly-nonlinear  (\v\  A>  1)  auto-oscillator  in 
the  intermediate  temperature  interval.  For  larger  temperatures 
or  smaller  oscillation  powers  po  (i.e.,  in  the  near-threshold 
region)  expansion  (19)  is  not  valid  anymore,  and  one  has  to 
analyze  the  full  equations  (1)  or  (5),  which  can  be  done  only 
using  numerical  methods. 

It  is  not  clear,  a  priori ,  which  of  the  above  derived  approxi¬ 
mate  expressions  for  the  ’’nonlinear”  auto-oscillator  linewidth 
(22)  or  (24)  is  a  better  approximation  for  the  description  of 
the  experimentally  measured  linewidth  magnitudes  in  strongly 
nonlinear  spin-torque  oscillators.  On  one  hand,  many  of  the 
experiments,  where  the  STO  linewidth  was  measured,  were 
done  at  a  room  temperature  T  ~  300K,  which  is  higher 
than  the  upper  limit  of  applicability  for  the  low-temperature 
expression  (22).  On  the  other  hand,  in  many  cases  the  ex¬ 
perimentally  observed  linewidths  of  spin-torque  oscillators  are 
rather  narrow  Aio  <  Tq,  so  that  the  conditions  of  applicability 
of  the  other  asymptotic  expression  (24)  are,  also,  violated.  As 
it  will  be  demonstrated  below,  the  direct  comparison  of  the 
currently  available  experimental  data  with  the  above  devel¬ 
oped  theory,  in  general,  demonstrates  that  the  low-temperature 
expression  (22)  gives  a  much  better  account  of  the  existing 
experimental  data.  However,  additional  measurements  of  the 
generation  linewidth  of  STO  performed  in  a  wide  range  of 
temperatures  are  necessary  to  clarify  this  important  question 
and  to,  possibly,  see  in  the  experiment  the  cross-over  from  the 
linear  to  a  squate  root  dependence  on  temperature,  when  the 
temperature  is  increased. 


D.  Near-Threshold  Region 

The  region  near  the  generation  threshold  r+(0)  ~  r_(0) 
(or,  for  a  spin-torque  oscillator,  I  ~  7th)  is  the  most  difficult 
region  for  the  theoretical  analysis.  In  this  region  there  are  no 
obvious  linearization  schemes  by  which  (1)  can  be  simplified, 
and  one  needs  to  employ  various  numerical  methods  to  evalu¬ 
ate  the  auto-oscillator  linewidth  near  the  generation  threshold. 

The  distortions  of  the  power  spectrum  S(tt)  in  the  vicinity 
of  the  generation  threshold  have  been  studied  in  [16]  using 
analytical  and  numerical  methods  based  on  a  non-stationary 
Fokker-Planck  equation  (5)  for  the  probability  distribution 
function  V(t,p,(j>).  The  results  of  this  analysis  show  that, 
in  the  near-threshold  region,  the  power  spectrum  S(Q)  does 
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not  have  a  simple  Lorentzian  form  (13a),  but  rather  can  be 
represented  as  a  sum  of  partial  Lorentzians 


sm  =  E 

3 


2A  UjPj 

(ft  -  uij)2  +  A  to2 


For  a  ’’linear”  auto-oscillator  (u>(p)  =  ujq  =  const)  the 
central  frequencies  u>j  for  all  the  partial  Lorentzians  are  equal 
uij  =  ojq,  and  the  overall  power  spectrum  S(Q)  remains 
symmetric  with  respect  to  the  central  frequency  ojq. 

In  contrast,  for  a  ’’nonlinear”  auto-oscillator  (oj(p)  ^  const) 
all  the  frequencies  ujj  are  distinct  and,  in  general,  distributed 
non-symmetrically  around  the  zero-temperature  frequency 
wg  =  oj{po).  As  a  result,  the  power  spectrum  S(Q)  becomes 
asymmetric  in  the  near-threshold  region.  Such  asymmetry  of 
the  power  spectrum  has  been  recently  observed  in  experiments 
with  spin-torque  oscillators  based  on  metallic  GMR  nano¬ 
pillar  [7], 

Another  specific  feature  of  the  behavior  of  a  nonlinear 
auto-oscillator  linewidth  in  the  near-threshold  region  is  an 
apparent  broadening  of  this  linewidth.  Like  in  the  case  of  a 
usual  inhomogeneous  linewidth  broadening,  the  sum  of  sev¬ 
eral  partial  Lorentzians  with  different  central  frequencies  ZJj 
results  in  the  power  spectrum  S(£l)  that  has  the  overall  width 
that  is  larger,  than  the  partial  linewidths  A ujj  of  individual 
Lorentzians.  This  effect  is  illustrated  by  Fig.  5a,  where  the 
linewidth  Aw,  obtained  from  a  single-Lorentzian  fit  to  the 
power  spectrum  S(Q),  is  shown  as  a  function  of  the  bias 
current  I  for  a  a  spin-torque  oscillator.  Good  agreement  is 
found  between  the  numerically  and  analytically  calculated 
linewidths  far  below  and  far  above  the  generation  threshold 
in  all  cases,  as  indicated  by  the  dashed  lines  representing  the 
linewidth  values  calculated  using  equations  (15)  and  (22).  In 
general,  the  overall  linewidth  A u>  decreases  with  the  increasing 
current,  except  in  the  region  near  the  threshold,  where  a  local 
maximum  of  the  linewidth  can  appear  in  the  case  of  a  strong 
frequency  nonlinearity  \u\  1.  Depending  on  the  magnitude 

of  the  nonlinear  frequency  shift  |jV|,  this  maximum  can  be 
several  times  larger  than  the  generation  linewidth  in  both 
below-  and  above-threshold  regimes  [16].  The  experimental 
measurements  of  the  power  spectra  of  current-induced  mi¬ 
crowave  oscillations  in  STO,  also,  show  the  line  narrowing 
followed  by  line  broadening  near  the  generation  threshold  with 
the  increase  of  the  bias  current  (see  Fig.  3c  in  [7]). 

Fig.  5b  shows  the  variation  of  the  generation  frequency  W, 
obtained  from  the  single-Lorentzian  fit  to  the  power  spectrum 
iS(fl),  with  the  bias  current  I  [16].  This  variation  is  relatively 
slow  below  the  generation  threshold  u;  «  ujq  and  exhibits 
a  fast  quasi-linear  change  with  current  above  the  generation 
threshold  due  to  the  nonlinear  frequency  shift,  in  agreement 
with  (1.30). 


E.  Comparison  with  Experiments 

The  above  presented  theory  of  the  generation  linewidth  of 
nonlinear  auto-oscillators  successfully  explains  a  number  of 
qualitative  features,  experimentally  observed  in  the  linewidth 
measurements  on  spin-torque  oscillators.  First  of  all,  the 
inverse  dependence  of  the  generation  linewidth  (22)  on  the 


Fig.  5.  Generation  linewidth  A w  (panel  (a))  and  frequency  shift  Su>  =  U—ujq 
(panel  (b))  as  a  function  of  the  bias  current  I  obtained  from  single-Lorentzian 
fits  to  spectra  computed  for  different  values  of  the  frequency  nonlinearity  N  = 
N/Fq  (from  [16]).  Noise  power  rj  =  1 0— 3,  nonlinear  damping  coefficient 
Q  =  2.  Dashed  lines  in  (a)  show  the  below-threshold  (15)  and  above-threshold 
(22)  linewidth  asymptotics.  Used  values  of  the  parameters  are  typical  for  in¬ 
plane  magnetized  spin-torque  oscillators  (note  that  the  nonlinear  frequency 
shift  N  is  negative). 


energy  of  oscillations  £(po)  explains  the  general  linewidth 
narrowing  with  the  increase  in  the  bias  current  and  oscillation 
amplitude  observed  in,  e.g..  Fig.  4  in  [17].  Note,  also,  that  the 
energy  £(po)  is  proportional  to  the  effective  volume  14ff  of  the 
auto-oscillator,  which  is  much  larger  for  a  spin-torque  auto¬ 
oscillator,  based  on  a  magnetic  nano-contact,  than  for  a  similar 
device  based  on  a  nano-pillar.  Thus,  the  linewidth  expression 
(22)  provides  a  natural  explanation  for  a  well-known  experi¬ 
mental  fact  (see  [18]-[20])  that  the  auto-oscillation  linewidths 
associated  with  devices  based  on  magnetic  nano-pillars  are,  in 
general,  several  times  broader  that  those  in  the  devices  based 
on  magnetic  nano-contacts. 

In  Fig.  6  we  show  the  comparison  of  the  generation 
linewidth  calculated  using  two  asymptotic  expressions  (22) 
and  (24)  with  the  results  of  experimental  measurements  of 
the  temperature  dependence  of  the  linewidth  of  a  spin-torque 
oscillator  performed  on  the  nano-pillar  devices  #1  (Fig.  6a) 
and  #2  (Fig.  6b)  in  [15]  (see  Fig.  2  in  [15]).  The  geometrical 
parameters  of  the  nano-pillar  devices  were  taken  from  [15]  and 
it  was  assumed  that  the  excited  magnetization  oscillation  mode 
is  pinned  at  the  pillar  lateral  boundaries  (see  [21]  for  details), 
the  Gilbert  damping  parameter  qg  =  0.01,  the  nonlinear 
damping  parameter  Q  =  3,  and  the  polarization  efficiency 
s  =  0.4  were  assumed  to  be  the  same  for  both  devices  [10]. 
As  it  is  clear  from  Fig.  6,  the  simple  analytical  expression  (22) 
obtained  in  the  low-temperature  limit  gives  a  reasonably  good 
estimate  of  the  observed  linewidths  at  different  temperatures 
for  both  nano-pillar  devices  with  the  same  parameters.  At  the 
same  time,  the  second  asymptotic  expression  (24)  obtained 
in  the  relatively  high-temperature  limit  significantly  overesti¬ 
mates  the  generation  linewidth.  Assuming  that  the  parameters 
of  the  two  devices  are  slightly  different  (which  is  possible  due 
to  different  nano-patterning  and  different  thicknesses  of  the 
’’free”  magnetic  layer),  one  can  obtain  much  better  quantitative 
agreement  with  the  experiment  [15]  using  the  low-temperature 
expression  (22).  Although  it  is  difficult  to  determine  from 
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Fig.  6.  Generation  linewidth  2Au>  of  a  spin-torque  oscillator  calculated 
from  the  approximate  expressions  (22)  (solid  lines)  and  (24)  (dashed  lines)  in 
comparison  with  the  temperature  dependence  of  the  linewidth  in  a  nanopillar 
device  #1  (a)  and  device  #2  (b)  measured  in  [15]  (black  dots).  Panel  (a)  shows 
the  linewidth  measured  at  the  second  harmonic  of  the  signal  2AuJ2  =  8Au). 


the  experimental  data  presented  in  Fig.  6  what  is  the  real 
temperature  dependence  of  the  generation  linewidth,  it  is  clear 
that  the  low-temperature  asymptotic  expression  (22)  gives  a 
much  better  estimate  of  the  generation  linewidth  magnitude. 

In  Fig.  7  we  compared  the  theoretical  (22)  and  experimental 
[17]  (see  Fig.  6  in  [17])  dependences  of  the  full  generation 
linewidth  2Acj  on  the  out-of-plane  magnetization  angle  Of,  of 
nano-contact  spin-torque  oscillator.  All  the  parameters  of  the 
nano-contact  device  Fig.  7  were  taken  from  [17],  while  the 
current  1  =  9  mA  and  magnetic  field  H0  =  9  kOe  correspond 
to  the  center  of  the  experimentally  studied  region,  and  the 
nonlinearity  parameter  of  positive  damping  was  again  chosen 
to  be  equal  to  Q  =  3  [10].  It  is  clear  from  Fig.  7,  that  the 
linewidth  dependence  on  the  bias  field  orientation  calculated 
using  the  low-temperature  asymptotic  expression  (22)  is  in 
good  quantitative  agreement  with  the  experimental  results 
from  [17],  In  contrast,  the  classical  result  for  a  ’’linear”  auto¬ 
oscillator  (18)  (see  dashed  line  in  Fig.  7,  which  shows  classical 
result  (18),  multiplied  by  10)  predicts  a  much  narrower  lines 
and  a  monotonous  decrease  of  the  linewidth  as  a  function 
of  the  out-of-plane  magnetic  field  angle  Of, .  The  ’’nonlinear” 
result  (22)  gives  a  reasonable  qualitative  and  quantitative 
description  of  the  experimentally  observed  behavior  of  the 
auto-oscillator  linewidth,  and,  in  particular,  it  predicts  the 
linewidth  minimum  around  9q  =  80°,  which  is  clearly  seen  in 
the  experimental  data.  This  linewidth  minimum  occurs  at  the 
’’linear”  magnetization  angle  0o,iin*  for  which  the  nonlinear 
frequency  shift  coefficient  N  vanishes  (see  Fig.  1. 6a),  and, 
therefore,  at  this  magnetization  angle  the  spin-torque  oscillator 
behaves  as  a  classical  ’’linear”  oscillator  with  the  generation 
linewidth  determined  by  (18). 

For  an  STO  based  on  an  anisotropic  nano-pillar  one  may, 
also,  expect  a  non-trivial  dependence  of  the  full  generation 
linewidth  2A u>  on  the  in-plane  magnetization  angle  (f>0,  since 
in  this  case  the  nonlinear  frequency  shift  coefficient  N  sig¬ 
nificantly  varies  with  <j> o  (see  Fig.  I.6b).  Recent  experiments 


Fig.  7.  Generation  linewidth  2A u  of  a  spin-torque  auto-oscillator  calculated 
from  (22)  (solid  line )  in  comparison  with  the  angular  dependence  of  the  nano¬ 
contact  STO  linewidth  measured  at  a  room  temperature  in  [17]  (black  dots) 
(from  [10]).  The  dashed  line  represents  the  classical  result  for  a  ’’linear"  auto¬ 
oscillator  linewidth  calculated  from  (18)  and  multiplied  by  10  . 

[22]  performed  using  a  spin-torque  oscillator  based  on  a  nano¬ 
pillar  geometry  where  the  ’’free”  layer  was  magnetized  by  the 
bias  magnetic  field  oriented  in  the  ’’free”  layer  plane,  indeed, 
demonstrated  strong  variation  of  the  oscillation  linewidth  with 
the  variation  of  the  in-plane  magnetization  angle  cf> o.  This 
behavior  of  the  generation  linewidth  finds  a  natural  explanation 
in  the  framework  of  the  asymptotic  linewidth  expression  (22). 
Fig.  8  shows  the  comparison  of  the  experimental  data  from 
[22]  with  the  angular  dependence  of  the  linewidth  2A u> 
calculated  from  the  asymptotic  equation  (22)  using  the  general 
expression  for  the  nonlinear  frequency  shift  coefficient  in  an 
anysotropic  magnetic  film  and  taking  all  the  other  parameters 
of  the  spin-torque  oscillators  from  [22].  It  is  clear  from  Fig.  8, 
that  the  linewidth  corresponding  to  the  magnetization  along 
the  hard  in-plane  axis  (<f>o  =  90°)  of  the  nano-pillar  ’’free” 
layer  is  much  smaller  than  the  linewidth  in  the  case  of  the 
magnetization  along  the  easy  axis  (cf>o  =  0).  It  is,  also,  clear, 
that  the  linewidth  expression  (22)  gives  a  good  qualitative  and 
a  reasonable  quantitative  description  of  the  experimental  data 
for  both  nano-oscillator  devices  used  in  the  experiments  [22]. 

We  would,  also,  like  to  note  that,  similar  to  the  case  of  an 
spin-torque  oscillator  with  an  out-of-plane  magnetized  ’’free” 
layer,  in  the  case  of  an  STO  based  on  an  anisotropic  in¬ 
plane  magnetized  nano-pillar  it  is  possible  to  choose  the 
direction  and  magnitude  of  the  bias  magnetic  field  in  such  a 
way  that  the  nonlinear  frequency  shift  coefficient  N  vanishes 
(see  Fig.  1. 6b).  Under  these  conditions  the  linewidth  of  the 
anisotropic  STO  will  have  a  minimum  possible  value,  corre¬ 
sponding  to  the  case  of  a  ’’linear”  auto-oscillator  (18). 
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